Probing freeze-out conditions in 
heavy ion collisions with moments of 

T 1 

^j: charge fluctuations 



3 

in 



Oh' 



F. Karsch 1,2 and K. Redlich 3,4 



1 Fakultat fiir Physik, Universitat Bielefeld, 



qj Postfach 100 131, D-33501 Bielefeld, Germany 

d : 2 

43 



Physics Department, Brookhaven National Laboratory, Upton, NY 31973, USA 

3 Institute of Theoretical Physics University of Wroclaw, 

PL-50204 Wroclaw, Poland 

4 Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
OO 

(N ; July 15, 2010 

o 

Abstract 
We calculate the first four moments of baryon number, electric charge and 
strangeness fluctuations within the hadron resonance gas model. Different mo- 
ments and their ratios as well as skewness and kurtosis are evaluated on the phe- 
nomenologically determined freeze-out curve in the temperature, baryon chemical 
potential plane. The model results and its predictions as well as relations between 
different moments are compared with the first data on net proton fluctuations 
in Au-Au collisions obtained at RHIC by the STAR Collaboration. We find good 
agreement between the model calculations and experimental results. We also point 
out that higher order moments should be more sensitive to critical behavior and 
will also distinguish hadron resonance gas model calculations from results obtained 
from lattice QCD. 

CERN-PH-TH-2010-161 



1 Introduction 

One of the major goals in theoretical and experimental studies of the thermody- 
namics of strongly interacting matter is to explore the QCD phase diagram at 
non-zero temperature (T) and non-zero baryon chemical potential (hb) DP- A 
central target of interest is the 'critical point' - a second order phase transition 
point, that has been postulated to exist in the T-/ib phase diagram [2j [3]. Al- 
though its existence is not yet firmly established, its presence would result in large 
correlation lengths, i.e. large fluctuations in various thermodynamic quantities. 
Remnants of these large fluctuations may become accessible in heavy ion colli- 
sions through an event-by-event analysis of fluctuations [I] in various channels of 
hadron quantum numbers (charges), eg. baryon number (B), electric charge (Q) 
and strangeness (S). In fact, at vanishing baryon chemical potential it has been 
shown that moments of charge fluctuations are sensitive indicators for the occur- 
rence of a transition from hadronic to quark-gluon matter [5]. They drastically 
change during the transition. 

When comparing theoretical equilibrium calculations of charge fluctuations 
with experimental results from heavy ion collisions, the crucial question is, whether 
at the time of hadronization the thermal system generated in these collisions has 
kept memory of the plasma phase or the expansion period during which it may have 
passed by a critical point [6]. If so, this may lead to an enhancement of fluctuations 
over ordinary thermal effects at freeze-out. On the other hand, it is well known that 
basic features of hadronization in heavy ion collisions are well described in terms of 
the hadron resonance gas model (HRG) [8]. The analysis of experimental data 
on the production cross sections of various hadrons in heavy ion collisions shows 
astonishingly good agreement with corresponding thermal abundances calculated 
in a HRG model at appropriately chosen temperature and chemical potential [9J 
[TO] . Abundances of strange and non-strange mesons and baryons are consistently 
described by freeze-out temperatures and baryon chemical potentials that are a 
function of collision energy only. 

If some memory of large correlation lengths persists in the thermal medium at 
time of freeze-out this must be reflected in higher moments of charge distributions. 
These moments are more sensitive to large correlation lengths and thus relax more 
slowly to their true equilibrium values at the time of freeze-out [UJ . It thus is an 
interesting question to test whether also more detailed information on the thermal 
distribution of hadron species, predicted by the HRG model, can be confirmed 
experimentally. In fact, one may consider HRG model results on moments of 
charge fluctuations as a theoretical baseline prediction for the currently ongoing 
low energy heavy ion runs at RHIC JT2] and future studies of charge fluctuations 
at the LHC. Any deviations from HRG model predictions would constitute evi- 



dence for new phenomena at the time of hadronization that have not been seen in 
experiments so far. 

Of course, eventually thermodynamics at freeze-out should be described by 
thermal QCD, eg. through lattice calculations. A direct comparison of exper- 
imental and HRG model results on higher moments of charge fluctuations with 
lattice calculations in the hadronic phase is possible [13], but is still difficult as 
so far most lattice calculations are performed with staggered fermions on rather 
coarse lattices. They need to be performed closer to the continuum limit to repro- 
duce the correct hadron spectrum [14} [T5] . Fortunately the calculation of ratios of 
moments is less sensitive to such cut-off effects [U[l3][l6]. At present, it seems that 
lattice QCD calculations of ratios of moments of charge fluctuations, performed at 
non- vanishing chemical potential by using a Taylor expansion of thermodynamic 
observables [171 [TH] , are in good agreement with HRG model calculations for tem- 
peratures below the transition temperature. We will consider this issue in more 
detail in Section 5 and indicate when such agreement may break down. 

In the following we will discuss the dependence of higher moments of fluctu- 
ations of baryon number, electric charge and strangeness on the collision energy. 
In particular, we will calculate ratios of quartic and quadratic (kurtosis), cubic 
and quadratic (skewness) as well as quadratic charge fluctuations normalized to 
their mean value along a phenomenologically determined freeze-out curve in heavy 
ion collisions. We also discuss correlations of different charges, eg. the correlation 
of baryon number and electric charge normalized to the squared baryon number 
fluctuations. 

In the next section we will summarize basic results on the parametrization of 
the freeze-out curve in heavy ion collisions. In section 3 we introduce moments of 
charge fluctuations and discuss their calculation in the HRG model. In section 4 
we compare the HRG model results with recent data on net proton number fluctua- 
tions in heavy ion collisions at RHIC energies obtained by the STAR Collaboration 
|19| and discuss additional fluctuation observables that may be analyzed to further 
characterize thermal conditions at freeze-out. Section 5 is devoted to a compari- 
son of the HRG model and lattice QCD calculations. We give our conclusions in 
section 6. 



2 Freeze-out conditions in heavy ion collisions 

Abundances of strange and non-strange mesons and baryons produced in heavy 
ion collisions in a wide range of collision energies are consistently described by 
freeze-out temperatures and baryon chemical potentials that are a function of 
collision energy only. In fact, the freeze-out curve T{^,b) in the T-jjlb plane and 
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Table 1: Parametrization of chemical potentials nx along the freeze-out curve 
using the ansatz given in Eq. [2j 



the dependence of the baryon chemical potential on the center of mass energy in 
nucleus-nucleus collisions can be parametrized by simple functions [20] 



T(n B ) = a - 6/i| - cfi% 



where a = (0.166 ± 0.002) GeV, b 
GeV~ 3 and 
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(0.139 ± 0.016) GeV" 1 , c 
d 



(1) 
(0.053 ± 0.021) 

(2) 



1 + e^/SNN 

with d = (1.308±0.028) GeV and e = (0.273±0.008) GeV -1 . This parametrization 
agrees with the phenomenological freeze out condition of fixed energy per particle 
of about 1 GeV |2T]. 

The energy dependence of strange and electric charge chemical potentials are 
obtained from the HRG model by demanding strangeness neutrality and isospin 
asymmetry in the initial state of Au-Au collisions. They can be parameterized 
in the same way as it has been done for the baryon chemical potential in Eq. 
The corresponding fit parameters are given in Table Q] and are shown in the left 
hand part of Fig. [U The right hand part of this figure shows the variation of 
baryon number density, n# = {Nb)/V, and electric charge density, uq = (Nq)/V, 
on the freeze-out curve. The strangeness density, ns = (Ns)/V, vanishes due to 
the imposed strangeness neutrality condition. We note that the ratio of baryon to 
strangeness chemical potential on the freeze-out curve shows only a weak depen- 
dence on the collision energy, 
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which is in agreement with findings on freeze-out conditions at RHIC |22j . 

Current results from lattice calculations [IH [151 E3] indicate that at \xb — the 
transition temperature, within errors, coincides with the freeze-out temperature 
extracted at the highest RHIC energy. However, although firm statements from 
lattice QCD need to wait for calculations closer to the continuum limit, it seems 
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Figure 1: The energy dependence of temperature T and chemical poten- 
tials for baryon number /i#, electric charge \xq and strangeness /is along the 
chemical freeze-out curve in units of GeV (left). The right hand figure shows 
the variation of baryon number density ub, electric charge density uq and 
strangeness density n$ along the freeze-out curve. 



that for [i b > the curvature of the transition line T([j,b) for fis > is distinct 
from the freeze-out curve [23]. This may indicate that at larger [ib the QCD 
transition and the hadronic freeze-out in heavy ion collisions are separated. The 
latter seems to occur at lower temperatures. 

If this separation of transition and freeze-out curves is indeed confirmed and/or 
hadronic fluctuations at freeze-out loose all memory of the passage through a 
possible transition region, then not only particle yields but also charge fluctuations 
should be characterized by thermal freeze-out conditions in the hadronic phase and 
may be well described by the HRG model. 

In the next section we will work out ratios of moments of various conserved 
charges on the freeze-out curve using the HRG model. In particular, we will 
calculate ratios of quadratic, cubic and quartic fluctuations of baryon number. 



3 Fluctuations in the HRG model 

The basic quantity that describes thermodynamics at non-vanishing chemical po- 
tential is the pressure. In the grand-canonical ensemble with baryon number, 
strangeness and electric charge conservation it is obtained from the logarithm of 
the partition function as 
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(4) 



In the HRG model the partition function contains all relevant degrees of free- 
dom of the confined, strongly interacting matter and implicitly includes interac- 
tions that result in resonance formation. The logarithm of the partition function 
in the HRG is obtained as a sum over all stable hadrons and resonances and their 
anti-particles, which can be separated into contributions from baryons and mesons, 

lnZ(T,^B,/iQ,/i 5 ) = J2 lnZ t( T ^Q^s)+ J2 lnZ r( T ,VB,VQ,Vs) 

i £ mesons i £ baryons 

(5) 
The contribution of each particle species of mass rrii, spin degeneracy factor gi, 
carrying baryon number Bi, electric charge Qi and strangeness Si is expressed as 

ln Zf(T, V, ft) = ^gimj £ V ^ K 2 (k mi /T) exp(fcc^7/T) , (6) 

fc=l 

where q = (Bi, Qi, Si), ft = (/j-b, MQ> Ms) an d K 2 (x) is the modified Bessel function. 
The upper and lower signs are for bosons and fermions, respectively. The first term 
in Eq. [6] corresponds to the Boltzmann approximation. 

The partition function (Eq. [5|) together with Eq. [6] provides the basis for the 
description of the thermodynamics of a system composed of hadrons and reso- 
nances being in thermal and chemical equilibrium. In the following we will focus 
on HRG model predictions on fluctuations of conserved charges and their higher 
moments. 

Fluctuations of the net charge N q and its higher moments are obtained from 
derivatives of the logarithm of the partition function with respect to the corre- 
sponding chemical potential. These derivatives define generalized susceptibilities 
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The first derivative, Xq , is related to the mean value M q of the net charge N q , 

M q = (N q ) = VT\ q V , (8) 



i.e. it is the charge density in units of T 3 , Xq = n q /T 3 . The second derivative 
with respect to \i q jT gives the variance a q = ((5N q ) 2 ) as 



I = VT^, (9) 



a q • -"• A.q 



1 Note that we introduce here dimensionless susceptibilities. This notation differs from 
that used in Ref. [19] by a factor T™" 4 . 



with SNg = N q — (Nq). The third x q an d the fourth \ q order moments can be 
expressed in terms of 5N q as 

((5N q f) = VT\f , (10) 

<(5Ag 4 > - 3((5N q f) 2 = VT^K (11) 

We also introduce skewness (S q ) and kurtosis (/%), which are generally used to 
characterize the shape of statistical distributions, 

9 ._<W> ,_ ((6N q )*) 

Z q - -3 , K q - -j ,3 . l^lZj 

u q u q 

Using Eqs.[8]to[T2]we may relate mean values, fluctuations, skewness and kurtosis 
to the generalized susceptibilities introduced in Eq. [7J 

rfi (2) (3) (4) 

u g _ Xq q _ Xq 2 _ Xq /-,q\ 

M a ~ J» ' Q<Jq ~ v (2) ' ^O - v (2) • ii<:ij 

y Xg Xg Ag 

Finally we also consider correlations of charges, which can be obtained as mixed 
derivatives of the pressure with respect to chemical potentials for charge X and 
Y, respectively, 

(nm) ^+>(r,/2)/r 4 

XxY d(^ x /T) n d(^ Y /T) m ' l ' 

We will discuss the behavior of the ratios Xbs iXb an d Xbq IXb f° r fn = 1, 3 
on the freeze-out curve. 

The above relations can be used to analyze properties of a thermodynamic 
system in equilibrium. If the evolution of a heavy ion collision results in an equi- 
librated state that lost all its memory of the previous evolution and hadronizes 
according to a hadronic statistical model, e.g. the HRG model, this will be reflected 
not only in particle yields but also in various moments of charge fluctuations. As 
higher moments become increasingly sensitive to large correlation lengths, any de- 
viations from HRG model predictions at freeze-out could indicate that the system 
memorizes that it went through a region with large correlation lengths. 

In the following we compare HRG model predictions with first measurements 
of the kurtosis, skewness and variance of net proton multiplicity distributions at 
mid-rapidity for Au+Au collisions at ^snn =19.6, 62.4 and 200 GeV performed 
by the STAR collaboration [19]. 



4 Charge fluctuations on the hadronic freeze- 
out curve 

In the Boltzmann approximation, which is a suitable approximation in the param- 
eter range considered, the HRG model provides a simple result for the thermody- 
namic pressure, 

fl = ^E^Kffl^K/T) cosh[( J B i /i i? + Sifis + QiVQ)/T] . (15) 

i 

The contribution of anti-particles is explicitly included in Eq. [15] through the 
cosh[x]-term. Thus, the summation is to be taken only over stable hadrons and 
resonances. 

A specific dependence of the pressure given in Eq. [15] on chemical potentials 
implies definite predictions for ratios of different moments of charge fluctuations. 
In particular, from Eqs. [7] and [15] it is immediately clear, that ratios of fourth 
{Xb ) an d second (x B ) order moments of the baryon number fluctuation as well 

(3) 

as the ratio of third (xb ) anc ^ ^ rs ^ or der moments are unity, 

(4) (3) 

Xb_ = 1 Xb_ = 1 r 16 \ 

(2) ' (1) ' ^ U ^ 

Xb Xb 

irrespective of the values of chemical potentials and temperature. A direct conse- 
quence of the above HRG model results is that simple relations hold for the mean 
net baryon number, kurtosis and skewness, 

kb&b = 1 ) k b M b = S b ctb ■ (17) 

The first relation is straightforward from Eq. [T3] the second is valid since 

- m - x b m b _X B X b _X B X b X b _ « n ns x 

X B a X B X B Xb Xb Xb 
In the HRG model the skewness Sb&b can be explicitly obtained from 

s _ Hi g baryons djjmj/T) 2 K 2 (mj/T) Smh[(nB + SjfMS + QiH Q )/T} 

Ei e baryons di{mi / T) 2 K 2 {mi / T) cosh[(/x B + Sms + QiHQ)/T] ' 
Consequently, for fj,g = \xq = one gets 

Sbctb = tanh(// B /T). (20) 
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Figure 2: The ratio of quadratic fluctuations and mean net baryon num- 
ber (<j b /Mb), cubic to quadratic (Sb&b) and quartic to quadratic (kb^e) 
baryon number fluctuations calculated in the HRG model on the freeze-out 
curve and compared to results obtained by the STAR collaboration [19]. The 
dashed curves show the approximate tanh(yU£/T) result for k,bct b and Sbct, 
respectively. 



This simple result arises from the fact that in the HRG model only baryons with 
baryon number B = 1 contribute to the various moments. 

In heavy ion collisions the strangeness and electric charge chemical potentials 
are much smaller than hb (see Fig. 1). The above relation thus can be considered 
to be a good estimate of skewness at chemical decoupling. We will show in the 
following that corrections due to non-vanishing electric charge and strangeness 
chemical potentials are indeed small for baryon number fluctuations. 



4.1 Comparison of the HRG model results on baryon 
number fluctuations with RHIC data 

The relations for skewness and kurtosis summarized in Eqs. (|17p . (|19p and (I20|) 
are generic results, expected to hold if thermodynamics is governed by the HRG 
model. Knowing the energy dependence of thermal parameters along the freeze- 
out curve (Eqs. [T]and [2]) we can directly verify if these particular relations, deduced 
within the HRG model, are consistent with recent findings of the STAR collabora- 
tion, which measured moments of baryon number fluctuations through net-proton 
number fluctuations |19|. 



Fig. [2] shows a comparison of the energy dependence of quadratic fluctuations 
(ct b ) normalized to the net baryon number (Mb), skewness Sb&b and kurtosis 
k b°'b obtained in the HRG model at chemical freeze-out with the STAR data. 

Obviously, the HRG model provides a good description of properties of different 
moments of net proton number fluctuations measured at RHIC energies. The 
reason for considering ratios of charge fluctuations rather than absolute values for 
different moments was, of course, that one is independent of definitions of the 
interaction volume and also is less sensitive to experimental cuts and systematic 
errors. Moreover, some of these ratios have an interesting interpretation, like e.g. 
the ratio xh /Xb which directly characterizes the dominant degrees of freedom 
carrying baryon number [5]. 

In addition it is, of course, also of interest to understand whether the 
HRG model can quantitatively describe the energy dependence of the STAR 
data [Pbjl on the first four moments, i.e. the mean, variance, skewness and 
kurtosigj. 

In order to compare the HRG model calculations with the experimental 
results presented in [19] we note that this analysis only explored fluctuations 
in a limited phase space. In fact, the data on mean particle yields differ from 
previous results obtained by the STAR Collaboration [22]. From Ref. [BI| 
one gets: M p _ p ~ 1.75 ± 0.25 and M p _ p ~ 3.5 ± 0.4 in the central (0-5%) bin 
of Au-Au collisions at v /J/v r /v = 200 GeV and 62.4 GeV, respectively. These 
values should be compared with M p _ p ~ 8 ± 1.8 and M p _ p ~ 15.4 ± 2.1 
obtained at mid-rapidity at corresponding energies in [22]. These data differ 
by a common factor K ~ 0.22. Part of the difference may be attributed 
to the fact that net proton fluctuation data in Ref. [TH] were taken in the 
restricted transverse momentum range 0.4 < p? < 0.8 GeV. 

In the HRG model used by us the thermal phase-space of all particles 
is not restricted. Consequently, in order to compare predictions for differ- 
ent moments of net proton fluctuations with experimental data one needs to 
rescale its thermal phase space by the above mentioned factor K ~ 0.22. Ef- 
fectively, this corresponds to rescaling the volume parameter VT 3 appearing 
for instance in Eq. [HJ by the i^-factor, although its origin is not necessarily 
related with a smaller volume of the system at chemical freeze out. 

The change of volume with energy on the freeze out line is calculated 
by comparing data on dN n -/dy at mid rapidity for different <Js with HRG 



2 Of course, as we have already verified consistency of three ratios with the HRG model 
calculations, only the energy dependence of one of theses observable provides additional 
information. 
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Figure 3: The energy dependence of mean (M p _p) (top, left), square root of 
the variance (cr p _p) (top, right), skewness (S p _p) (bottom, left) and kurtosis 
(Kp_p) (bottom, right) of net proton number fluctuations. The points are 
the RHIC results obtained by the STAR collaboration [19J . The lines are 
calculated in the rescaled (see text) HRG model on the freeze out curve. 



11 



model resultsjj We then obtain, V = [dN 7T ~/dy] data /n 1 ^I lG [T 1 j2], where in the 
HRG model the negatively charged pion density n^3 G is calculated using the 
relation between i/i and the thermal parameters given in Eqs. [I] and [2j Our 
results on V(y/s) extracted in this way are consistent with those obtained 
recently in Ref. [TO] . 

Fig. [3] (top left) shows the energy dependence of the first moment (M p _p) 
of net proton number in the HRG model with a volume parameter, VT 3 , 
rescaled by the factor K ~ 0.22. One can see in this figure that the HRG 
model results are consistent with the data. 

Taking into account the results for various ratios of moments shown in 
Fig. [2] it immediately follows that the HRG model will also describe the 
energy dependence of other moments, i. e. variance, skewness and kurtosis. 
These are also shown in Fig. [3j 

The good agreement of HRG model calculations with RHIC data allows 
us to make predictions for different moments of charge fluctuations covering 
the energy range of the RHIC low energy scan and the lowest energy for 
heavy ion collision at the LHC. We summarize the HRG model results at 
different energies in Table |2J 

4.2 Electric charge and strangeness fluctuations 

More subtle dependencies on temperature and baryon number arise in the 
case of electric charge and strangeness fluctuations where multiple charged 
hadrons get larger weight in higher moments and where meson as well as 
baryon sectors contribute. This results in characteristic deviations of the 
kurtosis, more precisely kq<Jq = Xq iXq > from unity and also the skewness 
no longer is simply related to tanh(/iB/T). In the case of electric charge 
fluctuations we may separate contributions of different charge sectors to the 
partition function. For instance, for n even, we may then obtain for moments 
of electric charge fluctuations, 

Xq ] = y^ (in Z\Q\=i(T, fi B , Hq,Hs) + 2 n In Z\ Ql=2 (T, fi B , fiQ, Us)) , (21) 

where the logarithms of partition functions, In Z\q\, are obtained from Eq. [5] 
by restricting the sums over mesons and baryons to the relevant charge sec- 



3 For a compilation of heavy ion data on charged pion yields at mid-rapidity at different 
energies see e.g. Ref. [8]. 
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Figure 4: The ratio of moments of electric charge fluctuations on the freeze- 
out curve (left) and correlations between baryon number and electric charge 
fluctuations (right). 

tors. From this it is obvious that deviations of kqOq = Xo /Xo from unity 
only arises from contributions of baryons with electric charge 2. Similarly 
the odd moments are modified. On the freeze-out curve this leads to a char- 
acteristic dependence of ratios of moments on the collision energy that is 
shown in Fig. HI In the energy range relevant for current low-energy runs at 
RHIC [12] as well as at LHC one has kqctq — 1.8, which varies only little 

With y/SNN. 

In addition one may analyze correlations between baryon number and 
different moments of charge fluctuations. Some results are shown in the left 
hand part of Fig. HI 

For completeness we show in Fig. [5] fluctuations and correlations in the 
strangeness sector of the HRG model. In practice it may be more difficult to 
compare this with experimental results as it will be crucial that the analysis 
allows for strangeness fluctuations in a sub- volume and will not impose the 
constraint of vanishing strangeness on event-by-event basis. 

5 Lattice QCD and the Hadron Resonance 
Gas 

The above considerations suggest that results on quadratic, cubic and quar- 
tic fluctuations of baryon number are in good agreement with HRG model 
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Figure 5: The ratio of moments of strangeness fluctuations on the freeze- 
out curve (left) and correlations between baryon number and strangeness 
fluctuations (right). 



predictions on the freeze-out curve at least for collision energies ranging from 
y/SNN = 200, 62.4 down to 19.4 GeV. In terms of thermodynamic parame- 
ters the currently analyzed range of collision energies at RHIC covers a rather 
narrow temperature regime, 160 MeV < T < 165 MeV. Changes in different 
moments of charge fluctuations with collision energy are mainly due to the 
variation of the baryon chemical potential along the freeze-out curve which 
covers the range 23 MeV < \i B < 210 MeV. 

To the extent that the freeze-out temperature for Hb/T ~ seems to 
be close to the transition temperature from hadronic matter to the plasma 
phase, the experiments cover a regime where 0.96<T/T C <1.0. In this narrow 
temperature range and for quark chemical potentials that are significantly 
smaller than T, i.e. [iqjT = /is/3T<0.44, there exists a lot of information 
on charge fluctuations from lattice QCD calculations. Due to the rather 
small value of fJ> q /T already lowest order Taylor expansions of quark number 
susceptibilities provide good guidance for the behavior of susceptibilities on 
the freeze-out curve [T6| [T8| [25] . Ratios of susceptibilities have been shown 
to be consistent with HRG model calculations even close to the transition 
temperature. To this extent the results obtained at RHIC are not only con- 
sistent with HRG model calculations but are also consistent with lattice QCD 
results. 

It may be conceivable that hadron fluctuations will only be sensitive to 
thermodynamics at freeze-out. One thus may ask whether fluctuations of a 
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completely thermalized system at chemical freeze-out in heavy ion collisions, 
bare any knowledge about a nearby critical point in the QCD phase diagram. 
In the following we will argue, that this is indeed the case. 

In fact, we know for sure from lattice calculations that at vanishing baryon 
chemical potential and for temperatures close to the transition temperature 
the HRG model has to fail in describing thermal moments of sufficiently 
high order. This reflects the existence of a chiral phase transition in QCD at 
vanishing light quark mass. It is a consequence of chiral symmetry restoration 
that moments Xb w ^ diverge at T c (hb = 0) for n > 6. At non- vanishing 
quark mass this leads to an oscillatory behavior of moments close to T c ; e.g. 
Xb w ill vanish at the transition temperature (TSJ [2S]- A direct consequence 
of the analytic structure of the QCD partition function in the transition 
region is that Xb w ^ become negative. This has been observed in lattice 
calculations [26] as well as in chiral models [27]. In the HRG model, on the 
other hand, all moments of baryon number fluctuations are positive. In fact, 
the ratio Xb iXb w ^ be unity in the HRG model for the same reason as 
Kfldg = 1. In particular at jib — 0, where the freeze-out curve is expected 
to be closest to the QCD transition curve or may even coincide with it, 
one should find strong deviations from unity. Lattice calculations suggest 
that Xb iXb vanishes at the pseudo-critical temperatures and rapidly rises 
for temperatures below, but close, to the transition temperature [25]. An 
analysis of these high order moments at LHC energies clearly would have 
the best chance to reveal differences between HRG model and lattice QCD 
calculations and to find evidence for critical behavior already at [Ib/T ~ 0. 

Current results on Xb /Xb are s ^ m noisy [25] but suggest that also this 
quantity may be well described by the HRG model up to temperatures close 
to the transition temperature. Such an agreement will stop to hold closer to 
the transition temperature. However, a better determination of this quantity 
in the transition region is necessary to quantify this and to confront it with 
experimental data. Eventually, it may be necessary to evaluate also Xb IXb 
ratio in order to observe striking differences between HRG model calculations 
and lattice QCD results in the energy range currently covered by RHIC and 
LHC experiments. 

Of course, if the critical point exist in QCD and if the hadronic freeze-out 
occurs within the critical region, then already the second moment of baryon 
number and electric charge fluctuations should deviate from the HRG model 
result. The higher order cumulants should exhibit even stronger sensitivity 
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\/$NN 


x$/x$ 


xflxf 


xgVxS' 


xg'/xg' 


Xbq IXb 


Xbq IXb 


7.7 


1.01 


0.99 


4.18 


0.49 


0.34 


0.79 


11.5 


1.05 


0.95 


5.39 


0.39 


0.32 


0.79 


19.6 


1.23 


0.81 


7.95 


0.27 


0.30 


0.77 


39.0 


1.87 


0.53 


14.25 


0.15 


0.28 


0.75 


62.4 


2.75 


0.36 


21.97 


0.09 


0.28 


0.74 


200.0 


8.20 


0.12 


67.80 


0.03 


0.27 


0.74 


2760 


111.1 


0.09 


922.4 


0.02 


0.27 


0.74 



Table 2: Ratios of moments of baryon number and electric charge fluctuations 
as well as their correlations for several values of the collision energy (in units 
of GeV) covering the RHIC low energy run and LHC (last row) energies. 



M) 



K Q°Q 



M) /,,(2) 



Furthermore, one has Kgd^ = xh IXb = 1 in the entire energy range and 



Xo iXq varies from 1.85 at low energies to 1.75 at high energies. 



to critical fluctuations showing larger deviations from the model predictions. 
New insight into this is to be expected to come from the first low energy run 
at RHIC which has been completed this year. 

6 Conclusion 

We have discussed properties of net charge fluctuations in nuclear matter 
within the hadron resonance gas (HRG) model. We have focused on the 
behavior of fluctuations related to baryon number, strangeness and electric 
charge conservation. Based on the phenomenological relation between ther- 
mal parameters and collision energy in heavy ion collisions we have calculated 
ratios of different moments of these fluctuations along the chemical freeze- 
out curve where secondary hadrons exhibit thermal and chemical equilibrium. 
We have also considered the energy dependence of baryon-strangeness and 
baryon-charge correlations in the HRG model. Our calculations covered the 
energy range of s/s NN > 4 GeV where the description of fluctuations within 
the HRG model formulated in the grand canonical ensemble is adequate. 

Establishing generic properties of fluctuations in the HRG model we have 
compared its predictions with the recent data of the STAR Collaboration at 
RHIC on baryon number fluctuations, expressed by net proton fluctuations, 
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determined in Au-Au collisions at three different collision energies. We have 
shown, that the STAR data are consistent with the HRG model results. The 
change of measured fluctuations with collision energy can be accounted for by 
the variation of the baryon chemical potential along the freeze-out curve. A 
list of results from HRG model calculations at collision energies that include 
the published STAR results, the parameters of the 2010 low energy scan at 
RHIC and collision energies at the LHC are given in Table |5J 

We have argued that the apparent agreement of the HRG model with 
the STAR data does not necessarily mean that in heavy ion collisions at 
chemical freeze-out the system has lost entirely its memory of the expansion 
period during which it may have passed through a region of the QCD phase 
diagram where correlation lengths are large, as expected in the vicinity of a 
critical point. If freeze-out occurs close to or in the transition region between 
hadronic matter and quark-gluon plasma this will, even in the absence of a 
phase transition, show up in higher order moments. The first four moments 
of baryon number density fluctuations agree well between HRG model and 
lattice QCD calculations even close to T c . Deviations form the HRG model 
results, however, have been seen in higher order moments. We have argued 
that in order to observe critical fluctuations in heavy ion collisions one would 
possibly need to measure even higher order moments. 

Our results on energy dependence of different moments and their specific 
relations in the HRG model can be used to characterize 'ordinary' thermal 
properties of higher order moments of charge fluctuations and to set a baseline 
for the observation of critical fluctuations in nuclear matter created in heavy 
ion collisions. 
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